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Abstract -The Brownian motion of a hot nanoparticle is described by an effective Markov theory 
based on fluctuating hydrodynamics. Its predictions are scrutinized over a wide temperature range 
using large-scale molecular dynamics simulations of a hot nanoparticle in a Lennard- Jones fluid. 
The particle positions and momenta are found to be Boltzmann distributed according to distinct 
effective temperatures Thbm and Tk- For Thbm we derive a formally exact theoretical prediction 
and establish a generalised Einstein relation that links it to directly measurable quantities. 
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I Hot Brownian motion is the Brownian motion of 
'^nanoscopic colloidal particles that have an elevated tem- 
S perature compared to their solvent [l] . The phenomenon 
^ is ubiquitous in modern biophysical and nanotechnological 
I— ^work, where one often relies on nanoparticles exposed to 
laser light as tracers, anchors, or self-propelling nanoma- 
^ chines [2||8j. In a number of innovative applications the 
.^heating of the particle is moreover deliberately exploited 
fT^for detection, manipulation and surgery on the nanoscale 
I>(9}{13]. In general, heat and vorticity diffuse much faster 
^^than the colloidal particles, typical difFusivities being on 
^^the order of lO^^m^s"^ versus 10~^^ . . . 10~^° m^ s"-^, re- 
^""^ spectively. On the basis of this characteristic Brownian 
'""^ scale separation, one may treat the solvent temperature 
. .and velocity distributions as stationary radial fields T(r), 
^ u(r) around the instantaneous particle position. To a 
k>( good approximation, the hot Brownian motion of a single 
; I particle can thus be characterised as a stationary nonequi- 
ty librium process. It is then natural to seek an effective 
equilibrium description in terms of Markovian stochastic 
equations of motion with effective friction and tempera- 
ture parameters [T 14 . 

Such effective parameters are certainly valuable el- 
ements of any quantitative description or application, 
but it is also clear that they will be less universal and 
more context-sensitive than their conventional equilibrium 
counterparts. For example, one has to expect different 



effective parameters for different degrees of freedom of 
the colloidal particle, such as translational and rotational 
motion, and they differ not merely by simple geometric 
or kinematic factors |I5) . Also momentum and confor- 
mational degrees of freedom turn out to be governed by 
distinct effective temperature, as already independently 
pointed out by Barrat and coworkers [16j. The best one 
can hope for is thus a systematic and quantitative under- 
standing of the origin of the effective temperatures and 
transport coefficients pertaining to the observables most 
relevant in practice. This hope is justified by the observa- 
tion that Brownian motion is a mesoscopic phenomenon 
and as such allows some coarse-graining over microscopic 
degrees of freedom. In other words, the effective temper- 
atures and transport coefficients of the colloidal particle 
should emerge form the "middle world" 
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of stochas- 
tic thermodynamics [T8| and fluctuating hydrodynamics 



19 -21 rather than directly from a much more intricate 
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microscopic description. As a consequence, the effective 
parameters may still be expected to be reasonably insen- 
sitive to many of the usually elusive (and often accidental) 
microscopic details, such as the precise functional form of 
the atomic interactions. 

While this is true in principle, the appropriate meso- 
scopic approach is not always entirely obvious and 
straightforward, a priori. Therefore, it is valuable to have 
direct access to a comprehensive microscopic characteri- 
sation of some model system. A standard way to achieve 
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this is via molecular dynamics (MD) simulations, which, 
in contrast to real experiments, provide complete control 
over the microscopic conditions. Following the pioneering 
work by Alder and Wainwright, there have been extensive 
investigations of the microscopic basis of classical fluid dy- 
namics in general, and of Brownian motion and its trans- 
port properties in particular 22-25 . In the following. 



we report on MD simulations of a hot nanoparticle in a 
Lennard- Jones solvent. We find that its Brownian motion 
is characterised by a set of four distinct effective temper- 
atures, two for its rotational and translational configura- 
tional dynamics and two for the corresponding momen- 
tum or kinetic (k) degrees of freedom. In particular, we 
demonstrate the following statements that hold for both 
rotational and translational Brownian motion (the explicit 
demonstration for rotational motion is deferred to a com- 
panion paper [Ts]). 

1. The characterisation in terms of effective tempera- 
tures for the configurational and kinetic degrees of 
freedom, previously investigated for a free particle 
[ll[l4)[T6], still holds in presence of potential forces. 
It is largely insensitive to microscopic details {e.g., 
to the solubility of the nanoparticle), but the "ki- 
netic temperatures" are sensitive to the precise heat- 
ing mechanism. 

2. The "configurational temperature" Jhbm may be ex- 
pressed in terms of the effective friction coefficient 
Chbm and diffusion coefficient I?hbm of the particle 
via a generalised Einstein relation. 



bJhbm 



= D 



HBM Chbm 



(1) 



This allows Thbm to be inferred from directly measur- 
able quantities. Over a broad range of particle tem- 
peratures, eq. ([T]) is found to be in excellent agreement 
with a formally exact theoretical prediction, eq. ([3|), 
derived in the supplementary online materials [26| 7^ 

Theory The derivation of the effective temperature 
Thbm that characterises the position fluctuations of a hot 
Brownian particle parallels the contraction of the fluctuat- 
ing Stokes problem, as established for equilibrium Brown- 
ian motion ^21^. In our analytical calculations, the solvent 
is idealised as incompressible and the temperature and vis- 
cosity distributions are approximated by radial fields cen- 
tred at the instantaneous particle position, which is well 
justified for many experimental realisations. The crucial 
quantity in the calculation is the dissipation function 



(r) = 27/(r)r(r) : r(r) 



(2) 



which gives the energy dissipated locally by the solvent 
fiow in terms of the strain rate tensor f = (Vu + Vu'^)/2 
for a given solvent velocity field u(r). Assuming for the 
fluctuating part of the solvent stresses the standard equi- 
librium correlations, evaluated at the local temperature 



T(r), our formally exact calculation, as detailed in the 
supplementary online materials [26], gives 



Thbm 



/^r(r)<^(r)d3 



(3) 



Analytical evaluation of eq. ^ for constant viscosity and 
thermal conductivity — corresponding to the standard ve- 
locity and temperature profiles 4u(r) = -3[U + r(U • 
r)]R/r - [U - 3f (U ■ r)]R^/r^ and T{r) = To + ATR/r, 
respectively — yields the exact result, Thbm = To+5Ar/12. 
Here, AT denotes the temperature difference between the 
ambient temperature Tq and the solvent temperature at 
the surface of the particle. The temperature dependence 
of the viscosity and the thermal conductivity gives rise 
to non-universal terms of higher order in the temperature 
increment AT. Their explicit evaluation requires some 
additional effort. Namely, one has to actually solve the 
generalised Stokes problem 



Vp= 2V-77r = 2rV77 + 277V-r, V-u = 



(4) 



for a radially varying viscosity field ri{r) to explicitly de- 
termine the velocity field u(r), subject to the usual no- 
slip boundary condition u(i?) = U. Although real flu- 
ids, such as water, are slightly compressible, the assump- 
tion V- u = is still valid, even if V • (pu) ^ 0, as the 
density variation is assumed to move stationarily along 
with the Brownian particle. From u(r) one then obtains 
(/)(r) and thus Thbm, and the effective friction coefficient 
Chbm of the hot Brownian particle. To solve eq. Q, a 
differential shell method has been developed that yields 
exact numerical results 
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For most practical applica- 
tions, this complication can, however, be sidestepped by 
using an analytical procedure that produces quite accu- 
rate predictions for (f){r) and Chbm from eq. Q |1,14|. 
If the temperature dependence of the solvent viscosity 
can be parameterised by the phenomenological relation 
ry(T) = ?7oo exp[^/(T - Typ)] (as e. 5., for water), and the 
density and thermal conductivity of the solvent can be 
assumed to be constant, it yields 

Thbm « To + ^ AT + , ^"^'^"/^-^ , AT^ , (5) 
" 12 22(To + AT-Tvf) 

which improves eq. (13) of Rings et al. (2010) [l] and cor- 
responding estimates in Rings et al. (2011) [14] (numerical 
errors less than 2% for AT/Tg < 1). For a more detailed 
discussion, including a corresponding expression for the 
Lennard- Jones fiuid and an explicit formula for the effec- 
tive diffusivity, the reader is referred to the supplementary 
online materials [26] . 

MD simulations In the simulations, the solvent and 
the Brownian nanoparticle are not treated as continu- 
ous media but are themselves made of atoms interact- 
ing via the Lennard- Jones radial pair potential U{r) = 



4e[(cr/7 



vl2 



(cr/r)°] truncated at r = 2.5a. The atoms 



belonging to the nanoparticle are additionally bound to- 
gether by a FENE potential U{r) = -0.5Ki?o log[l - 
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Fig. 1: Snapshot of the simulation: the Lennard- Jones atoms 
in the nanoparticle and the solvent are color-coded in order to 
visualise the temperature gradient. 



(r/Ro)^] with k = 30e/cr^, Rq = 1.5cr. As usual, we mea- 
sure lengths, times, and energies in terms of the Lennard- 
^Jma^/e and e, respectively. For liquid 

119.8K, 
We note that the 



Jones units cr. 

Argon, they correspond to cr = 3.405 A, e//cB 
m = 0.03994 kg/mol and r ~ 2ps 
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critical temperature for the bulk Lennard- Jones fluid with 
a cutoff of 2.5 is Tc = 1.186 28 . In our simulations, the sol- 



vent and the nanoparticle comprise 107233 and 767 atoms, 
respectively, corresponding to a length L « 51 of the peri- 
odic simulation cell and a particle radius i? « 5 (see fig. [T]) 
for a screenshot). Concerning finite-size effects, which are 
mainly due to the long range hydrodynamic interactions 
between the periodic image particles, we refer the reader 
to the supplementary online materials [26| , where we also 
give some details concerning the home-grown code and its 
parallel processing on graphics cards (GPUs). 

In a typical simulation, the system was first equilibrated 
in the NPT ensemble at the prescribed temperature of 
T = 0.75 and pressure of P = 0.01 using a Nose-Hoover 
thermostat and barostat. During the subsequent heat- 
ing of the nanoparticle to the temperature Tp, realised 
by rescaling the velocities of the atoms belonging to the 
nanoparticle in each time step, the system evolved in the 
NPH ensemble, while the particles near the boundary of 
the simulation box were maintained at T = 0.75. At 
least four independent trajectories of 2 • 10^ time steps 
St = 0.005 (corresponding to a physical duration of 100 
ns) were computed for each nanoparticle temperature Tp. 
Coordinates and momenta were recorded after attaining 
stationary thermal conditions. Note that different physi- 
cal heating mechanisms for the colloidal particle, such as 
a heat source residing within the particle (which does not 
directly affect its center-of-mass velocity) versus an ex- 
ternal heat supply (which generally does), correspond to 
different numerical heating procedures. 



Fig. 2: Radial temperature profiles in the cool solvent around 
a hot nanoparticle of radius i? fs 5 for the particle temperatures 
Tp = 1.00 •, 1.25 ■, 1.50 O, 1.75 A, 2.00 ▼, 2.25 O. Dashed 
and solid lines are fits obtained from Fourier's law for a thermal 
conductivity A=constant (dashed), as assumed in theory, and 
A(r) oc 1/T, as inferred from the fit to independent simulations 
of an isothermal bulk fluid (inset), respectively. 



Due to the finite compressibility of the Lennard- Jones 
fluid, the simulation does not comply with the idealiza- 
tions made in our theoretical calculations (incompressible 
solvent, constant heat conductivity). Therefore, to ac- 
curately describe the temperature profile in the Lennard- 
Jones solvent (excluding the discontinuity due to the Kap- 
tiza resistance at the particle surface [29]), we need to 
account for the temperature dependence of the thermal 
conductivity. Separate simulations of an isothermal bulk 
fluid along the state-space curve depicted in fig. [3] (dis- 
cussed below) yield an inverse temperature dependence 
for the thermal conductivity, A oc 1/T, as demonstrated 
in the inset of fig. [2j The main figure compares the corre- 
sponding solution of the heat equation, 



T(r) = To (1 + at/To) ^ 



(6) 



to the measured temperature profiles around the hot 
nanoparticle. Here, Tg 0.75 denotes the nominal am- 
bient temperature at infinite distance r ^ oo from the 
particle. 

While T(r) can be obtained by averaging the kinetic 
energy of the solvent particles in the vicinity of r, it 
is more subtle to deduce the local viscosity r]{r) from 
the available microscopic data. Under isothermal con- 
ditions, the viscosity can be computed from the micro- 
scopic stress tensor axy via the Green-Kubo formula rj = 
V{kQT)~^ fi^ {axy{t)axy{0)) dt. However, the strong tem- 
perature gradient in our system precludes the evaluation 
of the correlation function over a sufficiently large volume 
to get reliable results. We therefore evaluated the isother- 
mal viscosity t]{T) of a homogeneous bulk Lennard- Jones 
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Fig. 3: Sample points of the radial temperature and density 
profiles T{r) and p{r) around a hot nanoparticle maintained 
at various temperatures Tp, in the T - p plane. The filled 
circles ( ) correspond to states for which the viscosity ri{T) 
of the bulk Lennard- Jones fluid was determined in separate 
simulations. The inset compares the latter(D) to the empirical 
formula from eq. Q ( — ) and alternative expressions proposed 



in the recent literature 30 ( — ) |3l] (- • -) 



fluid with the intention to translate it into rji^r) = '7[T(r)] 
using T(r) from eq. ([6]). We observed that the radial de- 
pendence of the density p around the heated nanoparticle, 
which is neglected in the analytical calculations, can be 
quite substantial for the Lennard- Jones fluid. Plotting 
sample points of the radial temperature and density pro- 
files T{r) and p{r) around a heated nanoparticle in the 
T - p plane for various particle temperatures produced 
the state curve delineated by the data points in fig. [3) 
In determining the isothermal bulk viscosity rjiT) of the 
Lennard- Jones fluid using to the Green-Kubo formula, we 
therefore took care to vary the barostat pressure such as 
to confine the measured bulk states to this curve. In the 
inset of fig. [sj we compare our data for ri{T) to a phe- 
nomenological two-parameter fit 



\n[i^{T)h^] = {AlTY 



(7) 



Other observables that can directly be obtained from 
the simulation are the effective steady-state friction Chbm 
and diffusivity Dhbm of the hot Brownian particle. While 
good estimates for Drbm are deduced from the nanopar- 
ticle trajectories, determining the friction is slightly more 
subtle 32 33 . We inferred Chbm from the decay of the 



momentum auto-correlation function of the nanoparticle 
using the Brownian limit 22 . In the simulations, the 
measured force F on the colloid does not correspond to 
the random force ^ entering the Green-Kubo formula for 
the friction coefficient. From a generalised Langevin de- 
scription, it can be shown that the correlations (F(O)F(t)) 
and ($(0)^(i)) become equal only in the limit of a diverg- 
ing reduced mass p ^ oo. The Brownian limit amounts 
to first taking the mass of the colloidal particle to infin- 



Fig. 4: Distribution of the position coordinates of a hot Brow- 
nian particle in an asymmetric, anharmonic potential well for 
two nanoparticle temperatures. The Boltzmann distribution 
fits the data for the effective temperature Thbm (solid line) 
but not for Tk (dashed line), as it should. The inset depicts the 
distributions of the position (filled symbols) and velocity (open 
symbols) coordinates in a harmonic well of strength if = 50 and 
Boltzmann distributions for the effective temperatures Thbm 
(solid line) and Tk (dashed line). 



ity and then taking the thermodynamic limit of infinitely 
many solvent particles. As a consequence, the momen- 
tum relaxation time of the colloid becomes large compared 
to the typical relaxation time of the random force auto- 
correlations, and its momentum autocorrelation function 
becomes Markovian 32 . It is then found to exhibit an 
exponential decay 



(P(i)P(0)) = (P'(0))< 



-(Chbm/m)* 



(8) 



with the decay time determined by the friction coefficient 
Chbm and the reduced mass p. In practice, the Brow- 
nian limit is realised following the constrained-dynamics 



approach 34 



Knowing Chbm and -Dhbm, we can determine the con- 
figurational effective temperature Thbm using the gener- 
alised Einstein relation eq. ([T]). The temperature for 
the kinetic degrees of freedom is extracted from the equal- 
time velocity autocorrelation of the nanoparticle by ap- 
plying the equipartition relation p{\J^) = iksT]^- Note, 
however, that for hot particles, which are far from equilib- 
rium, the pertinence of the equipartition temperature 
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is a priori unclear. We refer the reader to Joly et al 
for a discussion of this aspect. 

Results and discussion In the remainder of this pa- 
per, we present our simulation results for the effective tem- 
peratures TjjBM and of translational hot Brownian mo- 
tion and test our theoretical predictions. 

To support the first claim formulated in the introduc- 
tion, we performed simulations in the presence of exter- 
nal potential forces derived from the harmonic potential 
V(r) = Kr'^ 12 for different values of the stiffness parameter 
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conditions such as the particle solubility or the physical 
realisation of the heating mechanism (internal/external), 
while Tk is sensitive to the heating mechanism. 

In summary, more than a century after Langevin pub- 
lished his famous equation 35 , we can now describe the 



overdamped Brownian motion of a hot nanoparticle by 



CHBMr=-VV(r)+^(t) . 



(9) 



0.3 0.4 

Ar I2Tq 

Fig. 5; Effective temperatures of translational hot Brownian 
motion. Simulation: Thbm (•) from the generalised Einstein 
relation, eq. ([TJ; apparent equipartition temperature Tk (O) for 
the particle velocity; solvent temperature at the particle surface 
(dotted) Theory: Thbm according to eq. Q (solid line) and the 
previous thermodynamic estimate, eq. (22) of Ref. ^14 (dashed 
line), both evaluated (assuming an incompressible fluid) with 
the exact numerical differential shell method |14|. 



K. Generalising the corresponding definitions for the free 
particle, the effective temperatures arc inferred from the 
particle's mean kinetic energy (fceTk = m(U'^)/3) and mean 
square displacement (fcBlHBM = K{r^)/3) in the harmonic 
well. In the supplementary online materials [26| , they are 
compared to results for a free Brownian particle main- 
tained at the same temperature to demonstrate perfect 
agreement to within the statistical errors. In the inset of 
fig.|4]it is moreover verified that the coordinates and mo- 
menta are indeed Boltzmann distributed with the respec- 
tive effective temperatures. The main plot demonstrates 
for two exemplary particle temperatures Tp that this gen- 
eralises to an asymmetric and anharmonic potential well 
V(r) = Y.i=i vixi) with v{x) = Kx^j'i + hx. Afiogether, the 
good agreement in all cases confirms that the particle ve- 
locity and position are distributed according to effective 
Boltzmann distributions governed by the same effective 
temperatures and Thbm a-s found for the freely diffus- 
ing particle. 

To establish also our second claim, fig.[5]provides a com- 
parison of our theoretical prediction for Thbm from eq. ([S]) 
and the simulation results obtained by means of the gener- 
alised Einstein relation, eq. ([I]). The good agreement val- 
idates eq. ([T]) and eq. (|3| over a wide temperature range. 
Also note that, at strong heating, eq. ([s]) fits the data sig- 
nificantly better than the previous thermodynamic esti- 
mate, eq. (22) of Ref. 14 — which amounts to a weighted 
average of 1/T instead of T in eq. ([s]), and is therefore 
systematically too small. As detailed in the supplemen- 
tary online materials [26| , the configurational temperature 
Thbm is found to be insensitive to the precise microscopic 



The effective Gaussian thermal noise ^(t) is characterised 
by {i{t)) = and {mUt')) = '^TuBuCHBuSit - t'). All 
coefficients are explicitly known over a wide temperature 
range, and accurate analytical expressions for them have 
been derived. Most importantly, we have established the 
following fundamental properties of the effective tempera- 
ture Thbm given in eq. (|3]): (1) it governs the Boltzmann 
factor for the probability distribution of a hot nanoparticle 
in an external potential V(r); (2) it may be expressed in 
terms of the directly measurable effective diffusivity T)hbm 
and friction Chbm via the generalised Einstein relation, 
eq. ([T]); and (3) it takes distinct values for translational 



and rotational degrees of freedom 15 . The momenta of 



the particle were also found to be Boltzmann distributed 
with and without external potential forces, but according 
to yet another effective temperature T^ that also takes 
distinct values for translational and rotational degrees of 
freedom, and, in contrast to Thbm, is sensitive to the pre- 
cise heating mechanism. Altogether, our findings strongly 
support the expectation [l] that the effective temperature 
Thbm can be treated as a bona-fide temperature for the 
configurational Brownian motion of an individual hot par- 
ticle, and thus pave the way for manifold applications. It 
remains as an intriguing open question how far this con- 
venient description can be extended to account for finite 
particle densities, anisotropics, [2j and self-thermophoresis 
[36]. 

Supporting Information A PDF file with technical 
and supplementary information containing a detailed for- 
mulation of the contraction of the fluctuating hydrody- 
namics problem to a Markovian Langevin description of 
the Brownian motion of the nanoparticle, convenient ana- 
lytical approximations for Thbm and T'hbm, a discussion 
of how to obtain Chbm and how to deal with finite size ef- 
fects in the numerical simulations, information about the 
parallel processing of the simulation code, and various sup- 
porting plots is available free of charge via the Internet at 
http://arxiv.org/| .. 
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